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A dependence of a relaxation rate on the shape of the demagnetized region for a longitudinal evolution of total 
magnetization to its equilibrium value following the ultrafast demagnetization is demonstrated. This shape- 
dependence is caused by a motion of the wave front inside the demagnetized region. The contribution of the 
wave front for spherically symmetric shape of the demagnetized region is up to 3 times and for cylindrically 
symmetric shape up to 2 times stronger than for one dimensional demagnetized region. This effect can be 
observed after the demagnetization by a tightly focused femtosecond laser pulse. 
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1. Introduction 

Recently, in the works where for a tightly focused (up to 1 pm) femtosecond laser pulse it was possible 
to observe the real picture of the spatial-temporal evolution under the action of a solitary pulse, the 
authors reported on the development of strongly inhomogeneous states. Along with the propagation of 
spin waves far away from the light spot (l] 13, the authors observed the formation of domains (3) and 
ring structures localized near the spot O. Strongly inhomogeneous distributions of the magnetization 
appear in the framework of superdiffusive mechanism of action of the laser pulse on magnetic metals 
mo and the existence of these distributions is estabhshed for magnetic heterostructures (6]I3- Thus, 
the analysis of the picosecond, longitudinal evolution of inhomogeneous, nonequhibrium distributions 
of magnetization has become increasingly important. 

The analysis of such states can be performed using the Landau-Lifshitz equation with a relaxation 
term proposed by Bar’yakhtar IMIOl. also referred to as LLBar equations (nj. LLBar equations were 
derived using general principles (Onsager’s relations, symmetry of the exchange interaction), and they 
have no alternatives for the study of the highly non-linear, nonequilibrium, non-uniform states with 
strong reduction of magnetization mentioned above. LLBar equations are well suited for a description 
of non-uniform states, such as magnetic solitons (12] ES) and Bloch points (HI, and give the explanation 
of the reversal effects These equations provide an explanation |7( of recent experiments O on 

magnetization recovery in laser-pumped Ni-Ru-Fe heterostructures, where the importance of the nonlo¬ 
cal character of the magnetization recovery is established (3. 

A longitudinal relaxation of a total magnetization of ferromagnets following the ultrafast demagne¬ 
tization within one dimensional (Id) model is determined by two effects: a homogeneous evolution of 
a magnetization inside the demagnetized region and the motion of the wave front from the outside to 
within the demagnetized region (T3. For a sufficiently large diameter of the demagnetized region, a ho¬ 
mogeneous evolution dominates. However, for a smaller demagnetized region, a relative contribution of 
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the wave front increases and for tightly focused femtosecond laser pulse the motion of the wave front 
can significantly enhance the relaxation I17II18I . 

In this paper, the effect of the size and shape of the demagnetized region on the relaxation rate of the 
total magnetization is analyzed. In order to demonstrate this effect, we consider the nonlinear, longitu¬ 
dinal evolution of the total magnetization following the ultrafast demagnetization for two limiting cases, 
namely, for the demagnetized region of cylindrical (CS) and spherical (SS) symmetries. Then, we compare 
these results with the results found within Id model. 

A strongly non-equilibrium state created by the laser pulse plays the role of the Initial condition for 
LLBar equations (7). For our case, since the motion of the wave front is associated with the change in the 
total magnetization of the ferromagnet, and the size of a transition region between demagnetized and 
non-perturbed regions is much larger than nanometers, the exchange (nonlocal) term, which retains the 
total magnetization of a sample, can be disregarded (compare with (7l[T3) and the main contribution to 
the equations of motion is determined by relativistic processes. Since during the relaxation of magneti¬ 
zation toward an equilibrium value the effective field is parallel to the magnetization, only the modulus 
of the magnetization M - |M| enters the equation. Adopting the Landau model for the free energy and 
introducing dimensionless variables, the nonlinear diffusion equation (NDE) is derived for the evolution 
of M 171 [181 

- m{l-, ( 1 . 1 ) 

dr 

where m(^,T) = MIMq is a dimensionless magnetization, Mo Is an equilibrium value (temperature- 
dependent) of the magnetization of a bulk material, ^ is a dimensionless radius in cylindrical or spherical 
coordinates measured In units of ro = ^2Ax\\ (as we consider cylindrically and spherically symmetric 
problems, only ^ enters the problem), A is an inhomogeneous exchange constant, - dM/df/ is a longi¬ 
tudinal magnetic susceptibility of a material in the equilibrium state and at zero magnetic field and t Is a 
dimensionless time measured in units of fo = 2 xi| ZyArMs, 7 is the gyromagnetic ratio, Ar is a dimension¬ 
less relaxation constant of the relativistic nature, Ms is the saturation magnetization. Simple estimates 
for nickel show that the value ro is of the order of a lattice constant, the characteristic time fo is of the 
order of a few picoseconds and the characteristic velocity ro/fo ~ 0.1 nm/ps = 100 m/s, see for details 
ITSl . Note that the use of the continuum approximation for distributions with characteristic sizes of the 
order of the lattice constant does not lead to qualitative errors, which can be seen from comparison of the 
results obtained numerically for discrete models and for their continual counterpart, see, e.g., (T^IZO). It 
is worth noting that this kind of NDE was first studied in the pioneering work by Kolmogorov, Petrovsky 
and Piskunov (ZT) and by Eisher (22) and stationary diffusive front propagations into unstable state have 
been found for this problem. 

After the pulse action, the demagnetized region (spot) is formed in a sample with a characteristic 
size of the order of a diameter of the laser beam. The value of magnetization is reduced inside the spot, 
m - mo < \, and outside the spot the magnet is non-perturbed, m - We consider two cases, when 
the demagnetized region has a cylindrical and spherical symmetries. To model such a situation, NDE is 
numerically solved for the following initial conditions 


m{^,T - Q) - mQ + 


_l-tttQ_ 

1 -I- exp [-4(^ - RfAla] ’ 


( 1 . 2 ) 


where Rq is the radius of the demagnetized region, a parameter a = (1 - m^) (dm/d( 5 ^)“^|m=(i+mo )/2 de¬ 
scribes the characteristic width of the transition region in the initial conditions. In the region of the ac¬ 
tion of the laser pulse (f < Rq], the magnetization tends to mo and outside this region (f > Jfo) m{£,,T - 0) 
tends to its equilibrium value 1. 

for a numerical analysis we consider the limiting case, when the laser pulse is focused to a diffraction 
limited spot of the order of 1 pm (2). The following values of the parameters (In the dimensionless units) 
correspond to the following situation: a radius of the demagnetized region (laser pump) Rq - 1250, the 
width of transition region a - 300 and the radius of the sample R - 1850. We consider the evolution 
of the initial state for the following minimal values of magnetization in the demagnetized region: 
mo - 0.9, 0.5, 10“^, 10“^, 0. The values mo - 0.5, 0.9 can be realized for a weak intensity of the laser 
pulse and correspond to the effective temperature of the spin system lower than the Curie temperature 
Tc- The values mo - 10“^, 10“^, 0 can be realized for a high power of the laser pulse. 
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2. Cylindrically symmetric initial distribution of magnetization 

Figure[^(a) presents the time evolution of the relative change of magnetization per unit length of the 
sample Mj; (t) 

AMz (T) ^ Mz (T) - Mx (0) 

AMx My, (+oo) — My (0) 

calculated in the time domain from t = 0 to r = T = 15 (in the dimensionless units to), for the above 
chosen parameters and for the initial conditions l |1.2| . Figure[^(b) demonstrates an example of the cor¬ 
responding evolution of magnetization m (^, r) on the time for mo - 10“^. 

Figure (a) shows that the fastest regime of relaxation of the total magnetization to its equilibrium 
value is realized for mo - 0.9. The relaxation time of AMj; (t)IAMy, as for Id model (13, decays with a 
decrease of mo and takes its minimum value at mo - 0. However, the dependence of AMy {t]/AMy on 
time is not linear for mo = 0 in CS case. The analysis of the evolution of m (i^, t) [figure[^(b)] demonstrates 
that similarly to Id model (T3IH), the evolution of the total magnetization to its equilibrium value is 
determined by two effects: the homogeneous evolution of the magnetization and the motion of the wave 
front within the demagnetized region. 

Let us derive an approximate analytical formula for the evolution of My (t) for CS case, which takes 
into account both the motion of the wave front and the homogeneous evolution of the magnetization and 
estimate the contribution of the wave front compared to Id model. Following (l3, the evolution of the 
wave front can be cast in the form 

m K, t) = mo (t) -I- [1 - mo (t)] mf« - - V't, , (2.2) 

where mo (t) does not depend on the coordinate and describes the homogeneous evolution of magnetiza¬ 
tion within the demagnetized region. The second term in the right-hand side of corresponds to the 
wave front, which moves within the demagnetized region with an amplitude 1 - mo (t) and the velocity 
V. The function mf describes the shape of this front; the value of mf tends to 1 if one moves outside the 
demagnetized region from the wave front and mf tends to 0 within the demagnetized region. Both V and 
the shape of the wave front depend on the distance from the center of the demagnetized region How¬ 
ever, this problem can be simplified by the following observation: figure[^(b) shows that the contribution 
of the wave front is essential only at the initia l sta ge of the evolution. This means that: (i) one can neglect 
the term (l/^)dm/d^oc l/(a-0 « and transforms to Id NDE; (ii) we can suppose that within 
the demagnetized region m (^, t) == mo (r) « 1 and the terms of the order of mo (t) « 1 can be neglected 
in this approximation. Then, the shape mf coincides with the shape of the wave front at mo = 0 in Id 



m(f, r) 





Figure 1. (a) The time evolution of AMy Ml AMy, derived from the numerical solution of NDE for the 
above chosen numerical parameters, different values of mo and for CS initial conditions. Full lines rep¬ 
resent numerical calculations and dashed lines are an approximation with the analytical equation (2.4| . 
(b) A corresponding to figure (a) evolution of the magnetization m{(, t) for mo - 10“^ on a time do¬ 
main from T = 0 to T = r = 15, taken step-by step after an interval of time T/IO. Dashed lines show the 
contributions of the wave front and the homogeneous evolution of magnetization. The arrow indicates a 
direction of motion of the wave front. 
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NDE and the velocity of the wave front V formed from the initial conditions can be found by equating 
the asymptotic of dD within the demagnetized region and that of the wave front:iY=a/4 + 4/ a (23l . For 
the chosen numerical parameters V ^ 75.01. Neglecting dispersion in O) and integrating this equation, 
the dependence of the magnetization on time mo (t) within the demagnetized region can be present as 
follows (7||: 

mo 

moM^—= - . (2.3) 

W nig + (l - ml) exp (-2t) 

Thus, at mo (t) « 1 the homogeneous evolution of magnetization leads to a decrease of the amplitude of 
the wave front [1 - mo (t)], and the motion of the wave front leads a decrease of the characteristic radius 
of the demagnetized region, which is equal to (i?o -Vt). 

Integrating dD over the unit length of the cylinder, the estimating formula for the evolution of M-z (t) 
is derived 

Mz (t) =: Mz (0) + nUlmo (t) + tt [1 - mo (t)] (2i?o - Vt) Vt, (2.4) 

where mo (x) is defined by (2.3| . Equation \2A) takes into account both the homogeneous evolution of 
magnetization and the motion of the wave front. 

Figure[^(a) shows a good agreement between the analytical equation l |2.4^ and the numerical results 
even in the regime, at which the magnetization within the demagnetized region mo (x) is not small, when 
both the shape and the velocity of the wave front depend on mo (x). The reason here is that for this regime 
the contribution of the wave front becomes smaU. 

Let us estimate the relative contribution of the wave front and the homogeneous evolution of magne¬ 
tization for different mo for CS case. We present 

AMz (X) dM™™ (x) dM^®* (x) 

^ ^^-, (2.5) 

QT QT QT 

where (x) and (x) are contributions of the wave front and the homogeneous evolution, re¬ 

spectively [these contributions are schematically depicted in figure[^(b) by dashed lines]. Differentiating 
j2.4| with respect to time and comparing the result with j2.5| , one can see that 

^^wave dM®®* 

—^— =^27i[Ro-VT)[l-mo[T)]V, ^ ==7i{Ro-VT)^[mo[T)-mo{T)^]. (2.6) 

dx dx 

Here, [1 - mo (x)] is the amplitude of the wave front, [Rq - Vt) is the characteristic radius of the demag¬ 
netized region at the moment x. 

Thus, (x)/dx takes its maximum value at the initial stage of the evolution. However, further, 

the the relaxation rate caused by the wave front decreases with a decrease of both the amplitude and the 
length of the wave front. 

The relative contribution of the wave front can be defined as (xl/AMj; at x — -i-oo, where 

^wave(x) jjg found by integrating |2.6| . Table presents M™^{t)/AMz for different mo, for Id 
model (^, CS and SS cases. Table [T] demonstrates that the contribution of the wave front for CS is 
stronger than for Id model (almost two times stronger for not small mo). To understand this observa¬ 
tion one should note that the wave front gives the main contribution at the initial stage. Then, estimating 
djffwave/(^^^dx) at X = 0 (initial stage), we see that this value is two times larger in CS case than that 
found within Id model. With a decrease of mo, increases and takes its maximum value 1 for 

mo = 0. 


mo 

0.5 

10"^ 

10-4 

0 

M^/AMzdd) 

0.05 

0.24 

0.51 

1 

M™''^/AMj;(CS) 

0.09 

0.41 

0.75 

1 

M™''^/AMj;(SS) 

0.13 

0.53 

0.85 

1 


Table 1 . The relative contribution of the wave front for Id, CS and SS cases for the above chosen values of 
the parameters. 
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To estimate the value of magnetization in the demagnetized region (for CS case), for which 
the relaxation rates for the wave front and the homogeneous evolution become comparable, we equate 
d^wave/dT= dM^^Vdr. This is a nonlinear, transcendental equation. However, a simple estimating for¬ 
mula can be derived by noting that the wave front gives the main contribution at the initial stage, that is 
« 1. Then, neglecting the terms like {rrig^]^, we get 

m^l^2V/Ro. (2.7) 

Thus, if the initial value of magnetization in the demagnetized region for CS situation mo « m^q, the 
wave front gives a substantial contribution to the relaxation of the total magnetization. Note, that is 
two times larger than for Id model ml^ == V/Rq flTl . 

3. Spherically symmetric initial distribution of magnetization 

The evolution of the wave front for SS case can be also presented in the form (2.2| . Using the assump¬ 
tions similar to the CS case and integrating over the volume of the sample, the estimating formulae 
are derived 

An In V^T^\ 

Mj_ (T) =: MxlO) -t —Rlmo{T) + AnVT[l-mQ{T)\ li?o -i?oUT-i- I' (3.1) 

where mo (t) is also defined by (H- Mx (t) is the total magnetization of the sample and 

— --47r(i?o-UT)^[l-mo(T)]U, ——— (i?o - Ut)^ mo (t) - mo (t)^ . (3.2) 

dr dr 3 

These formulae weU describe the results of the numerical calculations. Table [T] demonstrates that the 
contribution of the wave front for SS case is stronger than for CS case (almost 1.5 times for not small mo). 
Estimating dM™''®/(AM 2 ;dT) at the initial stage, we see that dM™''®/(AM 2 ;dT) at t = 0 for SS case is 1.5 
times larger than in the CS case, thus giving the explanation of the numerical data. Comparing different 
contributions in |3.2| we get 

mfq ^ 3V/Ro. (3.3) 

So, m^q is 1.5 times larger than for the CS case. 

4. Conclusion 

The relaxation of highly nonequilibrium, nonuniform states created by the femtosecond laser pulse 
in the relativistic relaxation approximation is governed by two scenarios: the motion of the wave front 
within the demagnetized region and the homogeneous evolution of magnetization within the demagne¬ 
tized region. The relativistic relaxation approximation is valid for widths of the transition regions be¬ 
tween the demagnetized and non-perturbed regions exceeding a few dozens nanometers (compare with 
(3). The contribution of the wave front for the SS shape of the demagnetized region is up to 3 times and 
for CS case is up to 2 times stronger than for Id model (l3; thus, the motion of the wave front within the 
demagnetized region leads to the dependence of the relaxation rate on the shape of the demagnetized 
region. For a sufficiently large diameter of the demagnetized region (laser spot), the homogeneous evo¬ 
lution dominates and the evolution of the total magnetization is the same for Id model, CS and SS cases. 
However, with a decrease of the diameter of the demagnetized region and with an increase of the power 
of the laser impulse, the relative contribution of the wave front increases and for mo < (m^q, m^q) the 
wave front can significantly enhance the relaxation and, consequently, the shape dependent effects can 
also be observed. This situation could be realized after demagnetization by a tightly focused femtosecond 
laser impulse. 
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no3AOB>KH5i eBO/iioi^m HaMarHineHOCTi (|)epoMarHeTi/iKa nic/in 
HaAUJBi/iAKoro posMarHinyBaHHB: po/ib posMipy ra (|)opMi/i 
posMarHineHoY o6/iacTi 

I.O. f1CTpeMCbKI/IW 

KhiBCbKuti HaL(ioHa./ibHMpi yHiBepcMTer iweHi Tapaca lileBteHKa, 

Byn. Bonofli/iMupcbKa, 64, 01601 Kmib, yKpama 

npofleMOHcrpoBaHO 3a./ie>KHicTb lubmakoctI penaKcattii nosflOBXHbo! eBO./iK)i4M hobhoT HaMarHiteHOcri (fiepo- 
MatHernKa nic.nfl HaALUBUflKoro poawarHiHyBaHHti Bifl c|)opMM posMarHineHoT o6;iacTi. Ifti sanexHiCTb oSywoB- 
.fiena pyxoM XBnabOBOro (fipoHiy bmh 6 posMarHiteHOi oSnacri. Bhccok XBi/i/ibOBOro c|)poHTy C(()epi/i4H0- 
CMMeTpuHHOi cfiopMM posMatHiteHOi o6aacTi flo 3 pas ra fl/w Ltn/iiHflpMtHO-cuMeTpuHHoT obaacri ftp 2 pas 
6iabujnpi, Hix flat! OflHOBUMipHoro BunaflKy. IJeM etfieKT Moxe cnocTepirarncb nicaa posMatHinyBaHHa ci/iab- 
HO ccjjOKycoBaHMM aasepHUM iwnyabCOM. 

KniouoBi cjioBa: HaflLUBi/iAKe posMarHi^yeaHHP, cpepoMarHeruK, <peMTOceKyHfl,Hi/ii^ jiaaep 
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